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A Dirac type variant of the xp Model and the Riemann Zeros 
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We propose a Dirac type modification of the xp-model to a x a ■ p model on a semi-infinite 
cylinder. This model is inspired by recent work of Sierra et al on the rcp-model on the half-line. Our 
model realizes the Berry-Keating conjecture on the Riemann zeros. We indicate the connection of 
our model to that of gapped graphene with a supercritical Coulomb charge, which might provide a 
physical system for the study of the zeros of the Riemann Zeta function. 



There have been an increasing interest in physical mod- 
els that might shed light on the Riemann hypothesis re- 
garding the zeros of the Riemann Zeta functions. This 
program aims at a possible realization of a conjecture 
by Polya and Hilbert on the zeros of the Riemann Zeta 
function [H, 0- The Polya-Hilbcrt conjecture [2J states 
that the complex zeros of the Riemann Zeta function on 
the critical line are described by the eigenvalues of a self- 
adjoint operator. One of the original proposals for such 
an operator arose from the work of Berry 0, Q| , followed 
by Connes Q and Sierra [B-GjJ with their collaborators. 
The attempts towards construction of such an operator 
areguided by various properties that are ascribed to it 
[H; EI i including the requirements that it should exhibit a 
chaotic classical dynamics with isolated periodic trajec- 
tories and that it should break time reversal invariance. 

Guided by these requirements, Berry and Keating Q 
suggested that the operator Hq = xp is a good candi- 
date for the realization of the Polya- Hilbert conjecture. 
The Bcrry-Kcating operator xp leads to chaotic classical 
dynamics and breaks time reversal invariance. However, 
it does not by itself leads to closed trajectories, which 
requires various additional identifications in the classi- 
cal phase space. The Berry-Keating operator was also 
analyzed by A. Connes 

An interesting physical realization of the Berry- 
Keating operator was given by Sierra and Townsend J7[, 
who observed that the xp-operator appears as the Hamil- 
tonian of the quantum Hall effect when restricted to the 
lowest Landau level [12]. They further argued that from 
this perspective, the Berry-Keating Hamiltonian has a 
quantum description consistent with the Polya-Hilbert 
conjecture. The relationship of Riemann zeros to other 
quantum systems such as the inverted harmonic oscil- 
lator [l3[ and the Morse potential [l4| have also been 
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discussed in the literature. 

Recently, Sierra and Rodriguez-Laguna [l(| have pro- 
posed a new operator Hi = x (p + y ^ where the particle 

is restricted to the (regularized) half-line l x < x < oo and 
lp is a constant. The operator Hi satisfies the essential 
properties of the Berry-Keating model and in addition 
gives rise to closed classical trajectories. In the quan- 
tum theory, the appearance of the - term in Hi leads 
to a nonlocal boundary condition. With an appropri- 
ate choice of the kernel of -, Hi was shown to admit a 
1-paramcter family of self-adjoint extensions with a spec- 
trum that is bounded from below and is thus a legitimate 
quantum Hamiltonian operator. For certain choices of 
the self-adjoint extension parameter, this model seems 
to be consistent with the Polya-Hilbert conjecture. 

The relation of the Hamiltonian operators in 0-0, 0, 
[ToL [IlT ] to the Polya-Hilbert conjecture appears through 
the Riemann- van Mangoldt formula 



N(E) = (N(E))+S(E), 



(1) 



where N(E) is the number of complex zeros of the Rie- 
mann Zeta function with positive imaginary part less 
than E, (N(E)) is the smooth part of N(E) and S{E) 
denotes a fluctuation in the counting (also known as error 
function) given by 



S(E) 



iE 



0(log(E)). (2) 



The smooth part of Eq. (fTJ) can be written as 



(3) 



where the phase $(E) of the Riemann Zeta function on 
the critical line (also known as Riemann- Siegel function 
0) is given by 

0(J5)=argr(i + ^-itf]og7r. (4) 
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The Hamiltonian operators introduced in 0, EH Eli 
provide an estimate for the phase $(E) of the Riemann 
Zeta function, which can be used in the Riemann-van 
Mangoldt formula ([T]) to make inferences about the zeros 
of the Riemann Zeta function. 

In this Letter, following the spirit of the work of Sierra 
et al [HI EH , we propose a variant of the xp model on the 
half-line where the operator p is replaced by a suitable 
two dimensional Dirac type operator p. The aim is to 
obtain a proper self-adjoint xp-operator on the half-line. 

Let us consider the Dirac operator 



a - p= O x p x + CTyPy, 



(5) 



where p a = —id a , with a = x,y, and a x ,a y are Pauli 
matrices. Thus, 





d x + id v 



d x - id y 




(6) 



This operator acts on two-component column- vectors val- 
ued on a Hilbert space defined on a semi-infinite cylinder. 
The semi-infinite cylinder is described by < x < oo and 
the y-dircction is a circle of radius R. There are a U(l) 
worth of possible boundary conditions that guarantees 
conservation of probability. This statement is equivalent 
to say that there is a U(l) family of self-adjoint operators 
ft on the semi-infinite cylinder. 

The above considerations on the operator ft should be 
contrasted with the fact that there is no self-adjoint op- 
erator p on the half-line. This statement can be checked 
by computing the deficiency indexes n± and then use 
the criteria discovered by Von Neumann [ID, EH- The 
deficiency indexes n± are defined as the dimension of 
the space of square- integrable solutions of the equations 
pip± = ±iip± on the half- line. For the operator p on the 
half line, we have n + = 1 and n_ = 0. Thus, by von Neu- 
mann theorem, there is no self-adjoint operator p on the 
half-line. Equivalcntly, there is no boundary condition 
on half-line that guarantees conservation of probability. 

We now consider the Hamiltonian is H = xp. This 
is our proposed modified version of the xp-model. This 
Hamiltonian acts on a suitable domain of the Hilbert 
space of square-integrable functions on the cylinder. Let 
us consider the normal ordered operator 1 



H = \fx ft \fx. 



(7) 



Similar to the operator ft, there is a U{\) family of self- 
adjoint operators H on the semi-infinite cylinder. Indeed, 
H inherits the same possible self-adjoint domains from 
ft. It may also noted that a different choice of normal 
ordering in the definition of H, like for instance H = 



1 In the final stages of this work, we learned that a very similar 
model to this one was already suggested by M. Asorey, J. Esteve 
and G. Sierra on unpublished notes. Their motivations was the 



(l/2)(xft + ftx) does not change the conclusions of the 
following analysis. 

For a two-component column-vector 



<J'i 

02 



the eigenvalue equation Hty = E^> leads to 

yfx {d x - idy) y/xip2 = iEip 1: 
\fx (d x + idy) \fxi\)\ = iEi/j 2 - 



(8) 



(9) 
(10) 



This system of equation is equivalent to a second order 
differential equation for one of the components, let us say 
tpi(x,y). From the topology of the y-dircction, we may 
consider the ansatz 



^2{x,y) = 



1 f{x)_ ^n+Hj 



/2tt x 
1 



27T 



X 



(11) 

(12) 



with n G Z and a G [0, 1). The parameter a is associated 
with the phase of the quasi-periodic boundary conditions 
along y-direction. It parameterizes the family of self- 
adjoint domains of the Hamiltonian. In more physical 
terms, a is related to a flux passing through the circle 
obtained by the compactification of the y-dircction. If 
we now set 



Ru 



2(n + a) ' 



we obtain 



d 2 cp 
du 2 



1 

2u 



E 2 



<P = 0, 



(13) 



(14) 



which is exactly of the Whittaker's form [17[. It may be 
noted that the expression for tp2 may be obtained from 
that of ipi . 

In order to set up the allowed boundary conditions, 
we first regularize the semi-infinite line < u < oo to 
uo < u < oo, for a positive uq- At uq we consider the 
boundary condition 



ip(u ) = 0. 



(15) 



same as ours. 



Later we discuss the crucial limit uq — > 0. The analysis 
of boundary conditions presented below can be cast in 
a more formal language of self-adjoint domains and the 
von Neumann theorem [Tsl . [l6| . The boundary condition 
on the wave-functions for the singular attractive inverse- 
square potential follo wing fro m the self-adjoint extension 
has been discussed in [18|, [l9| . The close relation between 
the boundary conditions following from the introduction 
of a cut-off to that obtained from self-adjoint extension 
has also been discussed in the literature pjl El, Hfll • Here 
we prefer to report on this route of regularization and 
renormalization, since it is more familiar among physi- 
cists. 
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The general solution of Whittakcr equation with 
boundary condition (fT5j) is 



A 

-W 



Wi_ iE {u ) Wi >+iE (u) 
h +isM Wi iE (u) , 



(16) 



where W{x) is Whittaker function 

W k ,m{u) = 



% u m+ 5 M ( m - k + i 1 + 2m; u 



with M(a, 6; m) being Kummer's confluent hypergeomet- 
ric functions (171 ]. 

We would like to have a solution that is square- 
intcgrable. For that, we first note that as u — > oo [17j ]. 



M(a, b; u) m y u 

1 (6 — a) 

Therefore, as u — > oo, 

1 



£(6) 
r(a) 



u 



<p(u) -¥ Ae 2 —= 
Ju 



-Wi 



w. 



;("o) 



r(l + 2iE) 



EK ~ UJ T{+iE) 

r(i - 2i£;)" 



r(-i^) 



(17) 



Thus the requirement of square-integrability implies the 
condition 



Wr_ iE (u ) _ r(l - 2iE) F(+iE) 
Wi t+iE (uo) ~ T(l + 2iE)T(-iEY 



(18) 



In our setup, the quantity u is a cut-off, which we 
would ideally like to take it to zero. Due to quantum 
instability associated with a strongly attractive inverse 
square potential, we cannot really remove it. As uq — > 0, 
we have fijl 



W a ,±iE(Uo) M| 



r±iE 



(19) 



Using Eqs. (fTS]) and (flT)]) . together with some basic prop- 
erties of Gamma function [17J , we get 



(?) 



2iE 



£(i 



+ f)r(| + f) 



f)r(f 



2 / 



Finally, using Eqs. ((H) and ([2U]). we get 



= 2ii9(.E)-KElnir 



"(I 



r(! + l #: 



(20) 



(21) 



This is comparable (but not the same) with equation (21) 
of 0. 

We have thus far shown that the xf model together 
with the boundary condition (TT5)) leads to an algebraic 
equation for the spectrum given by (|2Tj) . On one side, this 
algebraic condition contains the phase of the Riemann- 
Zeta function ^ which is part of the Ricmann-van Man- 
goldt counting formula ([1]). On the other side, apart from 



the ratios of certain Gamma functions, it contains an un- 
determined parameter, namely the cut-off u . 

The appearance of the cut-off uo in (fT5]) deserves some 
comments. As mentioned earlier, the cutoff uq is related 
to the self-adjoint extension parameter of the strongly at- 
tractive inverse square potential on the half-line [H, [lj| . 
Similar to the discussions in [ll| [2(J of a conformal quan- 
tum mechanics, this cut-off implies a breaking of the scale 
symmetry by quantization. That means a scale anomaly 
emerges in this problem. Indeed, from a classical per- 
spective, the xp model, and also the xp, is scale invariant. 
Also from another perspective, the Hamiltonian associ- 
ated with Eq. (fTl| in the short distance limit is singular 
and scale invariant. But the boundary condition (p~5|) 
breaks this scale invariance. The analysis of such singu- 
lar behavior for equations like Eq. (|14p is given in [19| . 
The works (2l| - [23| also contain important discussions for 
the emergence of scale anomaly in a similar context. 

The restoration of the scale anomaly at the quan- 
tum level is relevant for the problem of Ricmann Zeta 
function. Indeed, that allows us to enforce some of 
the requirements of Berry and Keating at the quantum 
level. One mechanism for the restoration of an other- 
wise anomalous symmetry is by the use of an appropriate 
mixed state 2AJ . This restoration mechanism in respect 
to the present problem is under current investigation and 
it will be reported elsewhere . 

Another objective of this Letter is to present, following 
the spirit of the work of Sierra and Townsend [7} , possible 
physical realizations for the H = xp model. 

To that end, let us first consider the Calogero model, 
whose solutions can be classified using the degree k of 
a polynomial that appears in the analysis (25j . It was 
shown in [26[ that for k = 0, the Calogero model with a 
complex coupling has real eigenvalues. The correspond- 
ing differential equation (see Eq. (12) of [26[) has ex- 
actly the same form as Eq. (|Ti|) , although the parame- 
ters appearing in the equations in these two cases have 
different interpretations. The main similarity between 
the two lies in the fact that they both describe a sys- 
tem of strongly coupled inverse square interaction. Such 
a system was first discussed by Landau, who observed 
that the strongly attractive inverse square coupling leads 
to a quantum instability that was characterized by "fall 
to the centre" [27| • Subsequently, a more complete quan- 
tum treatment was given by Case [l9| who found that the 
spectrum is unbounded from below and is characterized 
by a one-parameter family of an undetermined constant. 
In technical terms, this unknown constant is nothing but 
the self- adjoint extension parameter, which encodes in 
itself the combined effects of the short distance physics. 
In this work, as well as in (26|, the same role is played 
by the cut-off uq. Hence, the essential physics described 
here is that of a Calogero type system with a strongly 
attractive inverse square interaction. 

This naturally leads to the question if there is indeed 
a physical system which realizes such a potential. If a 
physical system exists whose eigenvalue equation is ex- 
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actly governed by Eq. (p~4|) . then the spectrum of that 
physical system would provide a concrete realization of 
the Polya-Hilbert conjecture. While we do not have an 
exact answer to this question, but there is a graphene 
system whose equations are again very similar to what 
has been described above. 

It is now well known that dynamics of low energy ex- 
citations in graphene is governed by a Dirac equation 
[28| . For our purpose, we consider the case of gapped 
graphene [29|, |3(|. In the presence of an external su- 
percritical Coulomb impurity I30l-l33i , the effective radial 
equation for gapped graphene [34| has the same formal 
structure as that of Eq. (fl"4]) . although again with a dif- 
ferent interpretation of the parameters. In addition, the 
solutions of a gap ped graphene system with supercritical 
Coulomb charge [34| have the same form as Eq. (TT4")) 
described above. When the external Coulomb charge 
in graphene is strongly attractive or equivalently in the 
supercritical region, the system again exhibits a quan- 
tum instability and "fall to the centre" (3ll - [36j . whose 
qualitative behaviour is similar to what we find here. 
In a graphene system, such an instability manifests it- 
self through characteristic features in the local density of 
states (LDOS), which in principle can be measured us- 
ing scanning tunneling microscope (STM) [3l| - [33| . With 
suitable rcinterpretation of the system parameters, the 
STM measurements of the LDOS in gapped graphene 
with a supercritical Coulomb charge may yield valuable 
information regarding the spectrum of our model consid- 
ered here. 

It may also be mentioned that the same equation ob- 
tained in f26j , appears in the analysis of non-relativistic 
AdS/CFT correspondence |37|. Typically, such strongly 
attractive inverse square operators arise also in the anal- 
ysis of near-horizon conformal structure of black holes 
[38l |39| . While none of these problems are identical to 
what has been considered here, the structural similarity 
of the equations and the corresponding solutions perhaps 
indicates a deeper relationship between these systems. 



The systems described above share a common feature 
that they all correspond to a strongly attractive inverse 
square potential characterized by quantum instability. In 
such systems, it is possible to employ a renormalization 
group scheme to address the issue of the instability (ioj . 
The instability indicates an incomplete understanding 
of the short distance physics. Typically such a system 
can be analyzed by first introducing a cut-off and then 
studying the renormalization group flow of the relevant 
parameters of the system, given by the corresponding 
/3-function. It would be interesting to study the conse- 
quences of the RG flow for the system described here. 

In summary, in this Letter, we have proposed a Dirac 
type variant of the H = xp given in Eq. ([7]). We then 
analysed its spectrum and related it with the Riemann- 
van Mangoldt counting formula for the complex zeros of 
zeta function Q. An important observation is that in 
order to find this quantum spectrum the scale invariance 
of the classical xp-model is necessarily broken. Based 
on this fact, we believe that the restoration of the scale 
invariance at the quantum level will have crucial im- 
plications to the Hilbcrt-Polya conjecture. Finally we 
suggested a physical realization of the present model in 
terms of the dynamics of the low energy excitations of 
the graphene in the presence of an external supercritical 
Coulomb impurity. 
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